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ABSTRACT 

We present a new chemodynamical code - Ramses-CH - for use in simulating the 
self-consistent evolution of chemical and hydrodynamical properties of galaxies within 
a fully cosmological framework. We build upon the adaptive mesh refinement code 
Ramses, which includes a treatment of self-gravity, hydrodynamics, star formation, 
radiative cooling, and supernovae feedback, to trace the dominant isotopes of C, N, 
O, Ne, Mg, Si, and Fe. We include the contribution of Type la and II supernovae, 
in addition to low- and intermediate-mass asymptotic giant branch stars, relaxing 
the instantaneous recycling approximation. The new chemical evolution modules are 
highly flexible and portable, lending themselves to ready exploration of variations in 
the underpining stellar and nuclear physics. We apply Ramses-CH to the cosmological 
simulation of a typical galaxy, demonstrating the successful recovery of the basic 
empirical constraints regarding, [a/FeJ-pe/H] and Type la/II supernovae rates. 
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1 INTRODUCTION 

The determination of elemental abundance patterns is one 
of the primary diagnostics of galaxy formation, with nu- 
merous spatial and temporal trends between age, kinemat- 
■ ics, and chemistry guiding our insights into the under- 
5^ ■ pining physical processes. Observations of abundance ra- 
, tios corroborate our understanding of the nuclear physics 
governing a-element production, in that they are pro- 
duced on sho r ter timescales than iron - peak e lements (e.g. 
Carbon et al] Il987l: lEdvardsson et al.l 1 19931 : iReddv et al] 
20061 : iRami'rez et'al]|2007l ). as a consequence of the mass- 
dependent nuclear burning processes acting within the rel- 
evant progenitor stars. Galactic chemical evolution (CE) 
models are predicated upon a coupling of these elemen- 
tal production sites/timescales with phenomenological (yet, 
empirically constrained) parameterisations of star formation 
and gas inflows/outflows. The resulting predicted abundance 
patterns can be compared directly with observations in order 



to shed hght on the formation and evolution of the system 
under study. 



The formalism associated with 
approach to galactic CE (e .g. 



the semi-numerical 



Talbot fc Arnett 1971: 



IPagel & Patchett^'l975VTinslev"l98(]';'Matteucci & Francod 
1989: .Carig i 1994; Gibson ,1997; Chiappini et al. 1993; 
I Ramirez et al.l |2007|) is a powerful tool when applied to 
sub-grid CE treatments within fully hydrodynamical sim- 
ulations. The inclusion of CE schemes has been achieved 
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based upon smoothed particle hydrodynamics (SPH). Key 
lessons can be learned from an examination of the role that 
CE plays in the physics of the interstellar medium. This 
is manifest in the metallicity-dependent radiative cooling 
rates of plasmas and their impact on the efficiency of metal 
transport throughout the disk an d its conse quent impact on 
stellar chemo-dynamics ( Scannapieco et al.ll2005 '). This im- 
pact upon turbulence-driven metal transport can be prob- 
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lematic, in light of known issues concerning the abihty of 
conventional treatments of SPH to resolve the associated in- 
stabilities in certain regimes; such problems are ameliorated 
(though not entirely) by Eulerian approaches to fluid dy- 
nami cs, including adaptive mesh refinement (A MR) sche mes 
(e.g lO'Shea et al.l l2005l : lAgertz et all l2007l : LTaskcr ct al] 
l2008l ). 

In its simplest form, interparticle 'mixing' of SPH par- 
ticles does not occur, i.e. metal-rich and metal-poor gas par- 
ticles may co-exist near each other without sharing/mixing 
of their associated metals. The impact of this lack of mix- 
ing is readily apparent in a galaxy simulation's metallicity 
distribution function and age- metallicity relation, as w ell as 
the abundance ratio plane (e.g. IPilkington et al.|[20l3 ). The 
inclusion of turbulent mixing models within S PH remedies 
this lack of implicit diffusion (|Shen et al.ll201(il l. even if the 
associated diffusion coefficient is a necessary additional free 
parameter (albeit, informed by turbulence theory). 

With the intention of providing a complementary 
(AMR) approach to extant (SPH) chemodynamical and 
semi-numerical CE models, we present what is, to our knowl- 
edge, the first cosmological AMR code which implements 
a temporally-resolved feedback and CE prescription. Writ- 
ten as a patch to the gravitational cosmological N-body 
and hydrodynamical code Ramses, we now include the ef- 
fects of Type H supernovae (SNell), Type la supernovae 
(SNela), and low- to -intermediate-mass asymptotic giant 
branch (AGB) stars, both from an energetic and chemical 
perspective. Nucleosynthetic processes are accounted for as a 
function of progenitor mass and metallicity. The CE module 
is described in ij2] The mechanics of marrying of this mod- 
ule to the AMR code Ramses, leading to the self-consistent 
chemodynamical code Ramses-CH, is outlined in 331 Fi- 
nally, in ^ we present a demonstration of the ability of 
the code to reproduce basic observational constraints. We 
emphasise that this Letter is primarily a methodological de- 
scription for coupling CE with AMR, and that future papers 
in this series will explore the response of Ramses-CH to the 
various input parameters and assembly histories. 



2 CHEMICAL EVOLUTION MODEL 

The underlying CE model used to determine the relative 
rates of SNeII:SNeIa:AGB, and the associated chemical en- 
richment, for a stellar population governed by a given IMF, 
is generated prior to the simulation being run. The result- 
ing look-up tables provide the SNell, SNela, and isotopic 
return rates as a function of time for a range of Simple Stel- 
lar Population (SSP) metallicities. The code (provided as 
part of the Ramses-CH patch) is flexible, allowing the user 
to readily modify relevant stellar physics, via the importa- 
tion of different SNe and AGB yields, as well as the IMF. 

For this first work, we have adopted a fairly stan - 
dard/conservative CE model, employing a iKroupa (|200ll ') 
IMF with stellar mass limits of 0.1 and 100 M©, re- 
spectively. We also used a SNIa delaye d time distribu- 
tion fo rmali sm similar to that p resent ed in iKobavashi et al] 
l|200(]h and iKawata fc GibsonI (|2003l ). with the simphfica- 
tion that the IMF slope for both the primaries and secon- 
darie s wa s taken to be identical. Stellar lifetimes were taken 
from iKodama fc Arimotd l|l997h and are dependent upon 
both mass and metallicity. The yields of SNell progenitors 
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Figure 1. Mass of elements ejected by stars as a function 
of initial mass. Also shown are the abundances for a sin- 
gle SNela (horizontal arrows) for comparison (the position 
along the abscissa is arbitrary, chosen to avoid conflict with 
other data). The mass above which stars are considered to 
be SNell progenit ors is indicated at 8 M(;). Data for AGB 
stars are taken from[v an den Hoek fc Groeneweg cn ( 199 7), SNell 
from Woo slev fc WeaveiT ifigbsV and SNela from llwamoto et al.l 
11999 ). Points connected by solid lines denote the original data, 
those connected by dot-dash lines show adopted extrapolations 
to lower masses. Extrapolations are linear and scaled to the mass 
of the progenitor star. 



(11—40 M©), f or metallicities spanning Population III to 
solar, are from 'Wooslev fc WeaveJ (j 19951 ): for m>30 Mq, 
we adopt the yields associated with the Model B explosion 
energi es, after iTimmes et al.l (|l995l ) and iKawata fc GibsonI 
()2003l ). Massive stars in the range 8—40 M© are assumed to 
explode as SNell, where the yields in the range 8—11 M© 
are found by scaling the elemental mass fractions of the 
IIM0 stars with the progenitor mass. The same process 
is used to calculate the yields of stars down to O.IM© 
using the lowest mass (O.SMg) A GB star available in 
va n den Hoek fc GroenewegerJ ([1993). The yields adopted 
are illustrated in Fig. [T] 

As not ed earlier, nucleosynth etic yields for SNela were 
taken from llwamoto et alj (Il999l ). In addition to the time 
constraints imposed by the mass range of the secondaries in 
SNela (binary) progenitors, , Hachisu et al. (1993,) also sug- 
gest the use of a metallicity 'floor' which suppresses the 
formation of low-metallicity ([Fe/H]< — 1.1) SNIa progeni- 
tors. In light of the ongoing controversy regarding this pu- 
tative metallicity floor, we have adopted the conventional 
assumption that low-metallicity binaries are capable of form- 
ing SNela progenitors. These SNela yields are also noted in 
Fig- H by the horizontal arrows. 

The time evolution of the isotopic ejection rate (per 
unit mass) from an SSP is shown in Fig.[2l The Ramses-CH 
chemical evolution model generates a family of such ejection 
rates for any combination of yield compilation and initial 
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Figure 2. Ejection rate of dominant el emental i sotop es, per unit 
stellar mass, as a function of age for a iKroupal ( 120011) IMF. The 
upper and lower panels correspond to solar and Population III 
metallicity simple stellar populations respectively. 

mass function. We have chosen to simply show the impact 
of the choice of one conservative combination of parameters. 
This should not be construed as implying that this combi- 
nation is necessarily the best or unique pairing; it is simply 
chosen to demonstrate the efficacy of the methodology. 



3 RAMSES-CH 

We have introduced the CE p rescription into the v3.07 pub- 
lic release version of Ramses (|Tevssieill2002l ). Prior to these 
enhancements, Ramses tracked the total gas metallicity (Z), 
under the assumption of the instantaneous recycling ap- 
proximation and treating Z as a passive scalar advected by 
the hydrodynamical flow. Our new chemodynamical version 
(Ramses-CH) also employs passive scalars in the tracking 
of the dominant isotopes of H, C, N, O, Mg, Ne, Si, and 
Fe, and the chemical composition of the gas from which 
the stellar particles f orm is 'tagged' onto the new parti- 
cles. As described by [Dubois fc Tevssieil (|2008l '). star par- 
ticles are created in the high-density gas {pgas > pth) and 
spawned by a random Poisson process following a rate given 
by p, — etpgas/tff,whevetff is the local free-fall time of the 
gas, (37r/32Gp5as)^''^, and e* the star formation efficiency. 
Each stellar particle enriches the surrounding interstellar 
medium (ISM) according to its individual chemical history 
recorded in the look-up tables described in ij2] depending 
upon the particle's initial mass and metallicity. Chemical en- 
richment and feedback processes are treated simultaneously 
through a kinetic feedback mode. While we could apply ki- 
netic feedback to all the stellar populations, we do so only to 
those stellar populations that include SNell events. A ther- 
mal feedback mode is used when the stellar population ages 
and enters an AGB and SNIa "phases" . The kinetic feedback 
mode aims at reproducing those expanding gas flows gener- 



ated by the collective explosions of massive stars. At each 
time-step, density, momentum, energy, and metals are de- 
posited into all gas cells situated within a feedback-'sphere' 
of a user-specified radius centered upon the young star par- 
ticles. We set up our SNell feedback-sphere radius to 2 grid 
cells. The velocity of the entrained gas, linearly interpolated 
with the radius, is given by the number of SNell events, 
the energy generated by each event tsNEsN {(-sn being the 
efficiency with which the energy Esn ~ 10""^ erg couples 
to the surrounding ISM), and the total amount of gas to 
be entrained by the bubble. The entrained gas includes the 
gas released/ejected during the SNell events as well as the 
gas swept up by the bubble (Dubois & Devriendt, private 
comm.). This latter component is parametrized as times 
the ejected gas mass [f-u, = 10 in this work, corresponding 
to a mass loading factor of r7,„ = 1 for a standard run with a 
massive star fraction of 10%). The galactic outflow scheme 
just described is not applied for the SNIa events and when 
the stellar population enters an AGB and/or SNIa phase, 
feedback processes and chemical enrichment are handled 'lo- 
cally' (thermally) and confined to the gas cell within which 
the stellar particle sits. The feedback algorithm is simplified 
by this approximation; future developments of the code will 
include kinetic SNIa feedback. 

This new chemodynamical version of Ramses was then 
applied in the generation of a multi-resolved, cosmologi- 
cal simulation of a late-type disk galaxy (henceforth called, 
109-CH) with a virial mass of 7.1 x lO^'^ M©, whose ini- 
tial conditions match those outlined in an earlier study 
(|Sanchez-Blazquez et al.ll2009l ). To remind the reader, the 
box size for this run was L=20 Mpc, with a dark mat- 
ter particle mass of 6 x 10^ Mq in the most refined region. 
The spatial resolution, kept roughly constant in physical 
size during the simulation, reaches —436 pc at z=0, 

with lmax=16 levels of refinement. For this resolution, we 
choose a star formation density threshold pth corresponding 
to 0.3 cm~^ and an efficiency of 1%. The same SNe energy 
coupling efficiency of esjv = 1 is assumed for both SNell and 
SNela. A polytropic equation of state T — Tth{p/ pth)"''^ is 
used in the high-density areas with a temperature threshold 
of 2900 K and a polytropic index 7=2, allowing the Jeans 
length to be resolved by more than 4 cells at all times. Note 
that our feedback parameter choices are suitable for the im- 
plementation described here at resolutions of 100-500pc. 



4 CHEMODYNAMICS OF AN L. GALAXY 

We introduce this new grid-based chemodynamical tool 
to the community. We now demonstrate the efficacy of 
Ramses-CH through a presentation of the chemical proper- 
ties of 109-CH. Future papers in this series will extend be- 
yond this initial demonstration into a comprehensive chem- 
ical tagging and chemodynamical exploration of a suite of 
higher-resolution cosmological disks spanning a range of en- 
vironments and assembly histories, realisations of which are 
described in iPilkington et all l|2012l) . 

The first metric to consider when introducing SNell 
and SNela in parallel with a relaxation of the instanta- 
neous recycling approximation is the predicted supernovae 
rates and a comparison with empirical constraints. We show 
the time evolution of these rates in Fig. [3] For 109-CH, a 
present-day SNela rate of 0.131 SNuM (SNe per century 
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Figure 3. The star formation rate for 109-CH is sliown in black 
(refer to tlie left-hand ordinate) and the corresponding SNell and 
SNela rates are shown (refer to the right-hand ordinate) in blue 
and red, respectively. Note that the SNII rate is not precisely 
proportional to the star formation rate as would be the case for 
data simulated using the "standard" version of RAMSES. The SNII 
rate at each time is now dependant on the star formation rate of 
past as well as present time-steps. 



Figure 4. Abundance ratios of stars in a disk region of galac- 
tocentric radius between 5 and 11 kpc and within 3 kpc of the 
disk plane. Particles are coloured according to their age. Obser- 
vational d ata is plotted in g rey, triangles are very metal-poor 
stars from lCavrel et aP 2004h. diamonds are thick disk and halo 
stars from Gratton ct al. I ll2003l '), plusses are disk dwarf stars from 
Edvardsson ct al. (199^. All data has been normalised to the so- 
lar abundance determination of lAnders fc Grevessd l ll989l ). 



per 1O^°M0 stellar mass) was found, and a SNell rate of 
0.959 SNuM. Both the absolute values and SNelLSNela ra- 
tio (^ 7) are consistent with those found bv lMannucci et all 
l|2008l ) for field Sbc/d galaxies: 0.140+o o3g SNuM for SNela 
and 0.652t°;i^* SNuM for SNell. 

Moving beyond the SNe rates, the abundance ratios of 
readily observed elements are regularly employed to con- 
strain the timescales of star formation, and therefore both 
feedback and fundamental nucleosynthesis. The recovery of 
empirical trends found locally in the solar neighbourhood 
is a necessity for any CE model. Such observations demon- 
strate a clear correlation between a-element and iron abun- 
dances, in the sense of their being an a-enhanced plateau 
at lower metallicities (below, say, [Fe/H]'^— 0.7) with a sys- 
tematic decline to solar values, at higher me t allicities (e.g 
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200d ). This empirical behaviour in [a/Fe] — [Fe/H] can be 
seen in Fig.|4l where, in this case, we are showing the obser- 
vational trends for [O/Fe] vs [Fe/H]. Also shown in Fig. U is 
the same distribution of [O/Fe] — [Fe/H] for the star particles 
at redshift z=0 within an analogous 'solar neighbourhood' 
for the simulated disk 109-CH. We should emphasise that 
the solar normalisation emplo yed for both simul ated and 
observed data is that of i Anders fc Grevessd l|l98£j ). 

A rigorous analysis of how the choice of SNeH, SNela, 
and AGE yields, in addition to the IMF and SNela pro- 
genitor model impact upon the chemodynamical evolution 
is left to future papers in this series, but it should be 
clear that even with this first 'test', the qualitative chem- 
ical properties are not inconsistent with observations of the 



local plateau-|-decline behaviour seen in the [a/Fe] — [Fe/H] 
plane. Abundance ratios are recovered and the qualitative 
behaviour of the [Fe/H]~— 0.7 knee is also seen in the simu- 
lated data. The knee-feature seen in Fig. |4]at [O/Fe]=0.15, 
[Fe/H]=-0.2 is attributed to the bursty star formation pro- 
file which naturally creates multiple knee features as the 
SNn:SNIa rate fluctuates. 

Examination of the age distribution (represented by the 
colour-coding shown in the inset to Fig. lU reveals rapid 
early enrichment in [Fe/H], similar to the age-metallicity 
relations predicted by classical CE models. Specifically, it 
takes ~3 Gyr to reach a metallicity [Fe/H]~— 0.4, driven by 
the initial phases of intense star formation, after which the 
age-metallicity relation fiattens and the rate of growth of 
[Fe/H] consequently slows (even while SNela are becoming 
more important). This phase is characterised by the abun- 
dance ratio 'strata' see in Fig. |4l with discrete 'arcs' appear- 
ing at decreasing values of [O/Fe] as time progresses. 

The infiuence of the cosmological environment of this 
simulation is apparent in the abundance patterns of the 
galaxy. The sub-sample displayed in Fig. |4] exhibits the sig- 
nature of merger events, e.g. the feature at [O/Fe]=0.15, 
[Fe/H]=-0.2. The full extent of this is only apparent when 
examining all stars in the galaxy where discrete 'streams' 
with chemical properties distinct from the rest of the galaxy 
are seen. These arise from accretion of satellites that have 
lower [Fe/H] values and remain chemically distinct. Larger 
mergers may have a similar abundance to the primary galaxy 
but bring gas that can reignite a quiescent galaxy and ac- 
celerates the production of Fe in the short term. 
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5 SUMMARY 

We present a new chemodynamical simulation code that pro- 
duces feedback to account for long-lived stars and the ele- 
ments that they produce. The more sophisticated SN feed- 
back scheme improves the kinematic properties of the stellar 
fraction (this will be detailed in future work) and gives ac- 
cess to additional constraints on the sub-grid physics. It is 
clear that the galactic CE for our galaxy does not per- 
fectly reproduce Milky Way observations, however we be- 
lieve it serves to demonstrate the validity of this approach. 
Future work using this code may merit the inclusion of cut - 
ting edge nucleosynthesis models (e.g. iDohertv et "al]l2O10| ) 
and will make a comprehensive study of the influence of the 
parameters involved in the underlying CE model. 
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